Complete quantum control of the population transfer branching ratio between two 
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A five-level four-pulse phase-sensitive extended stimulated Raman adiabatic passage scheme is 
proposed to realize complete control of the population transfer branching ratio between two degen- 
erate target states. The control is achieved via a three-node null eigenstate that can be correlated 
with an arbitrary superposition of the target states. Our results suggest that complete suppression 
of the yield of one of two degenerate product states, and therefore absolute selectivity in photochem- 
istry, is achievable and predictable, even without studying the properties of the unwanted product 
state beforehand. 



I. INTRODUCTION 

Optical control of atomic and molecular processes has 
attracted great theoretical and experimental interest in 
recent years 0, 0. The seminal study of Brumer and 
Shapiro Q showed that the relative phase of two inde- 
pendent excitation pathways between the same initial 
and target states can be optically manipulated to gen- 
erate interference control over the ratio of populations 
tranferred to degenerate product states. However, this 
weak-field picture of the controlled quantum dynamics is 
not very useful in the case of strong fields, where the num- 
ber of interfering pathways that relate the initial state to 
the same target state can be enormous. 

The stimulated Raman adiabatic passage (STIRAP) 
£| method for population transfer is a powerful strong 
field approach to the control of atomic and molecular pro- 
cesses. In its simplest three-level version, STIRAP gen- 
erates coherent population transfer between two states 
of the same parity without populating the intermediate 
state. The strong field eigenstates are conveniently rep- 
resented as time-dependent superpositions of the three 
field-free states. The key element in the STIRAP pro- 
cess is the creation, via counter-intuitively ordered pump 
and Stokes pulses, of a null eigenstate that has a node on 
the intermediate field-free state and correlates with both 
the field-free initial and target states. There have been a 
number of successful extensions of STIRAP to multi-level 
systems H I 0, 1 1 M M EMS M M US 13 111 
In particular, Kobrak and Rice [llj (KR) have proposed 
a five-level extended STIRAP scheme to achieve prod- 
uct selectivity in photochemistry. This scheme was later 
shown to be the essence of the important Chen-Shapiro- 
Brumer (CSB) strong field approach [H H E(j to the 
control of photodissociation reactions (note that the CSB 
approach was experimentally realized in the Na2 system 
|20|). Since then, there has been a renewed interest in 
studies of controlled molecular dynamics from a STIRAP 

perspective |21l l22l l23l 124 l25l I26L m nn m mm 

The central interest in the KR scheme is to control the 
population transfer branching ratio between two degen- 



erate target states. This issue is of great importance to 
selective photochemistry because photodissociation reac- 
tions typically form degenerate product states, with each 
product state correlated with a different potential energy 
surface. Unlike the degeneracy associated with, for exam- 
ple, different magnetic sublevels where the exact degen- 
eracy can be lifted by Zeeman splitting 9] and certain 
sublevels may be addressed individually by taking ad- 
vantage of the selection rules associated with polarized 
laser fields 0, 0] , the exact degeneracy between differ- 
ent chemical product states in the continuum cannot be 
lifted. Then, if one of the degenerate product states is 
coupled with a particular state by a laser field with an 
arbitrary polarization, so is the other product state. In- 
deed, neither the KR scheme nor the CSB strong field 
approach offers a means for completely suppressing pop- 
ulation transfer to an arbitrarily chosen product chan- 
nel. A simpler and potentially very useful approach to 
control of population transfer that takes advantage of a 
large lifetime difference between two degenerate product 
states has been proposed j2^|, but the control is incom- 
plete when both target states rapidly decay. 

Recently, we advocated a measurement-based ap- 
proach that may be used to completely shut off the yield 
of one of two degenerate product states [3(| • In particu- 
lar, we have shown that a strong measurement of the pop- 
ulation of the branch state in the KR scheme modifies the 
spectrum of the system such that the nonadiabatic cou- 
pling and the quantum interference between two partic- 
ular adiabatic states can be manipulated. This approach 
to population transfer control overlaps with the study of 
the quantum Zeno and anti-Zeno effects. However, this 
approach does not totally resolve the issue of achieving 
complete control of the population transfer branching ra- 
tio. In fact, the typical yield of the desired product state 
is not 100% and analytical conditions for the complete 
suppression of the yield of one product state remain to 
be established. 

In this paper, we propose a five-level four-pulse ex- 
tended STIRAP scheme that offers complete control of 
the population transfer branching ratio between two de- 
generate product states. The results strongly suggest 
that complete selectivity of product formation is in prin- 
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ciple achievable. As shown below, in addition to intro- 
ducing one more laser pulse than in the KR model and 
in our measurement-based approach, the price that is 
paid for complete control is that the relative phases be- 
tween the four laser pulses must be well maintained. This 
requirement is experimentally demanding, but with the 
evolution of laser technology we expect it to be met. Even 
if the relative phases between the four laser pulses cannot 
now be controlled, we expect that a number of attractive 
properties inherent to our five- level four-pulse scheme will 
motivate new experiments. For example, we demonstrate 
that one does not need to analyze the properties of one of 
the two degenerate product states in order to predict the 
complete suppression of its yield. We also demonstrate 
that an arbitrary superposition of the two degenerate tar- 
get states can be created. Finally, we show that complete 
suppression of the yield of one of two degenerate product 
states can be realized when the product states have life- 
times that are short compared with the duration of the 
laser pulses. 

The achievement of complete controllability in multi- 
level systems may be rigorously formulated as a math- 
ematical problem. Some existence proofs that establish 
sufficient conditions for complete control have been dis- 
covered jH El HI- The Huang-Tarn-Clark theorem Hjj 
is believed to the strongest result, but it applies only 
to systems with discrete and nondegenerate states. Ra- 
makrishna et al. showed that for a quantum system with 
a Hilbert space of dimension L, the necessary and also 
sufficient condition for complete controllability is that 
the field-free Hamiltonian and the interaction Hamilto- 
nian induced by a control field generate a Lie algebra of 
dimension L 2 j35j. This criterion is applicable to degen- 
erate multi-level systems, but the required computations 
to generate the commutators of the Lie algebra struc- 
ture can be demanding for large L, vary drastically from 
system to system, and provide no hint as to how a con- 
trol field can be constructed. Even if the existence of 
complete controllability in a multi-level system that in- 
volves degenerate states is established using a Lie alge- 
bra analysis, that result is unlikely to provide hints for 
the following: (i) how to design a conceptually simple 
control scheme to realize the complete control; (ii) how 
a physically valid approximation changes the picture of 
controllability, (iii) how the decay of product states af- 
fects the controllability, and (iv) how the controllability is 
influenced by constraints on the system. Consider, for ex- 
ample, the generation of a complete population transfer 
pathway in a five-level system that never populates more 
than two levels. This control scheme can be formulated 
as a mathematical question that is very different from the 
original question of how to achieve complete controllabil- 
ity in the same system. Indeed, as suggested by an early 
study by Shapiro and Brumer [37|], such a population 
transfer pathway is very unlikely to exist (but not neces- 
sarily impossible) because the dimension of the Hilbert 
subspace that is disallowed by the constraint is too large 
relative to that of the rest of the Hilbert space. Never- 



theless, as we show below in a straightforward manner, 
a complete population transfer pathway in a five-level 
system that does not populate three levels does exist, in- 
dicating that the achievement of complete controllability 
may be possible with reasonable approximations (such as 
the adiabatic approximation and the rotating-wave ap- 
proximation) even though the underlying mathematical 
question seems difficult to resolve or suggests the oppo- 
site. Hence, our five-level four-pulse control scheme is 
also of great theoretical interest in the context of the for- 
mal question of complete controllability without or with 
strong constraints. 

This paper is organized as follows. In Sec. II we 
describe, in connection with the original KR five-level 
model, our five-level four-pulse phase-sensitive extended 
STIRAP scheme for the complete control of the popu- 
lation transfer branching ratio between two degenerate 
product states. In Sec. Ill we show that our control 
scheme can also be used for the creation of arbitrary su- 
perpositions of two degenerate target states. Some sim- 
ple computational examples that support our theoretical 
analysis are presented in Sec. IV. Our results are dis- 
cussed and set in context in Sec. V. 



II. A FIVE-LEVEL FOUR-PULSE EXTENDED 
STIRAP SCHEME 

A. The model system 

A schematic diagram of our five- level four-pulse scheme 
is shown in Fig. ^ The initial state |1) is coupled with 
the intermediate state \2) by a pump pulse, state |2) is 
coupled with two degenerate product states 13) and |4) 
by a Stokes pulse. As in our previous work [ill IT^. we 
have adopted the common (but less rigorous) treatment 
wherein the product states |3) and |4), if embedded in the 
continuum, are described as discrete levels with complex 
energies whose imaginary parts give the decay rate con- 
stants. Unless specified otherwise below, we first assume 
that the lifetimes of all the five levels are much longer 
than the duration of the laser pulses. This assumption 
is not essential for the success of our control scheme, but 
is convenient in the description of the physical picture of 
the control mechanism. Due to the degeneracy of states 
1 3) and |4), alteration of their population transfer branch- 
ing ratio (denoted by B) cannot be achieved unless other 
laser-induced couplings are introduced. This observation 
motivated Kobrak and Rice to introduce a "branch state" 
1 5) that is also coupled with states |3) and |4) by a third 
"branching pulse" . Here we consider also a fourth laser 
pulse, called the "control pulse" , that further couples the 
branch state with the initial state. From Fig. Q]it is seen 
that the control pulse introduces an evident symmetry to 
the system: state |5) can be thought of as a second inter- 
mediate state, and the control (branching) pulse plays a 
similar role to that of the pump (Stokes) pulse. As shown 
below, this extension from the KR model to our five- 
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level four-pulse scheme leads to a new class of STIRAP- 
like processes that can be advantageously manipulated 
by laser phases. 
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FIG. 1: A schematic diagram of the five-level four-pulse 
phase-sensitive extended STIRAP method for the complete 
control of the population transfer branching ratio between 
two degenerate product states |3) and |4). The initial state 
is |1). Qp, Qc, fis3 and £2s4, f2s3 and ils4 represent the 
Rabi frequencies associated with the pump pulse, the control 
pulse, the Stokes pulse, and the branching pulse, respectively. 
Under certain conditions 100% population can be transferred 
from state |1) to state |4) without ever populating states |2), 
3), or 1 5). Cases (a) and (b) represent two possible level 
configurations. 

We assume that all the laser fields are on resonance 
with the respective transitions and that the pulse shapes 
are Gaussian. The electric fields of the pulses are given 
by 



E x {t) = E x f x (t) cos(cu x t + (fix), 
with the envelope function 
fx(t) = cxp 



<t-t x f/T x ] 



(1) 



(2) 



In Eqs. (JTJ and J5J, lo x is the frequency of the laser 
pulse, (fix is the associated phase, T x is the pulse width, 
tx represents the timing of the peak electric field am- 
plitude Ex, and X — P,S,B,C for the pump, Stokes, 
branching, and control pulses, respectively. Below we 
confine ourselves to level configuration (b) in Fig. ^ the 
other configuration can be analyzed in a similar fashion. 
Using the rotating wave approximation (for discussions 
on the validity of this approximation, see Ref. |23j) and 
the interaction representation, the system Hamiltonian 
is given by 



(3) 
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Here the six Rabi frequencies (see Fig. ^)) are 



rip = 


|/i 12 |(£ P /2)/p(t)exp[i( 


^p4 


"12)], 


(4) 


n S3 = 


| M23 |(S s /2)/ s (t)exp[i( 


<fis + 


"23)], 


(5) 


sls4 = 


\ m \(E s /2)f s (t)exp[i( 


<fis + 


"24)], 


(6) 


fi B 3 = 


\fi 35 \(E B /2)f B (t)exp{- 


■i((t>B 


- "35)], 


(7) 


Q B 4 — 


M(E B /2)f B {t) exp[- 


■i{<f>B 


- 045)]) 


(8) 


n c = 


|/i 15 |(S c /2)/c(t)exp[i( 


<fic 4 


"15)], 


(9) 



where and a k j represent the magnitude and the 

phase of the transition dipolc moment fi^j between states 
|fc) and \j). 



B. The null eigenvector 

It is quite unexpected and fortunate that the eigen- 
values and eigenvectors of H, denoted by A& and |Ajt) 
(k = 1 — 5) below, can be obtained analytically. In par- 
ticular, we find that H has one null eigenvalue (denoted 
Ai), with the associated eigenvector given by 



|Ai> 



1 



^S4^* B3 — ^S3^B4 







o* o* 



(10) 



where N\ is the normalization factor. Except for one 
accidental case discussed below, the other eigenvalues are 
found to be nonzero, with A 2 , (k = 2 — 5) satisfying a 
simple "quadratic" equation, namely 

= {xtf-x 2 k n 2 + |n ss nj 4 -fi^n^l 2 

n* c n Si \ 2 + \npn* B3 - n* c n S 3\ 2 (n) 



I n B ^ 



where il 2 is the sum of the absolute squares of all the six 
Rabi frequencies. Using Eq. the analytical forms 

for |Afe) (k — 2 — 5) can also be obtained. However, since 



4 



states |Afc) (k — 2 — 5) are in general not degenerate with 
|Ai), the nonadiabatic coupling between \\ k ) (k = 2 — 5) 
and | Ai ) can be made negligible by use of sufficiently 
strong fields. With the weak nonadiabaticity assump- 
tion, only the null eigenstate is relevant to theoretical 
consideration of the dynamics of population transfer. 

Consider first two special cases in which our scheme 
reduces to the original KR model where only three pulses 
are used. In the first case fic = 0, fip ^ 0, and state |5) 
is the branch state. Then one has 



|Ai) 



fis4fij3 3 — ^S3^*B4 


o* o* 
— o* o* 





(12) 



STIRAP scheme is a marriage between phase-insensitive 
STIRAP and coherent phase control. Indeed, comparing 
Eq. (0 with Eqs. {BJ and {T^Jl, one sees that the null 
eigenvector |Ai) in Eq. (|10|l may be understood as a co- 
herent superposition of the two null eigenvectors \0)kb. 
and \0)' KH associated with the KR model. That is, in 
essence, we are able to utilize two KR schemes simulta- 
neously by adding one more laser pulse to the system. 
Our five-level four-pulse approach is also related to and, 
as seen below, more ambitious than, a recent proposal by 
Karpati and Kis to achieve selective adiabatic population 
transfer to degenerate magnetic sublevels [l7|. 



C. Complete yield suppression of one product state 



The associated branching ratio is given by B — 
\tt B4 \ 2 /\n B3 \ 2 = |/i4 5 | 2 /l^35| 2 . Similarly, if fi P = 0, 
Oc 7^ 0, then state |2) plays the role of the branch state 
and 



|Ai) - \0)'kr = 





— rigClgi 




(13) 



The associated 

I2/IO 12 _ 



value of B is then given by 
tostf/Ps*? = Irf/lrf- Clearly, the value of 
B in either case is determined entirely by the properties 
of state 1 5) or state |2), is a constant during the popula- 
tion transfer, and does not carry any information about 
the laser phases. Further, for each desired value of B one 
has to seek a different branch state and then adjust the 
laser frequencies, or, alternatively, detune the branching 
pulse from the branch state in an effort to generate the 
desired value of B. Unless one of the transition dipole 
moments /135, /i45, /123 and /i24 happens to be zero, 
complete suppression of the yield of one product state 
or varying the value of B from to 00 is generally out 
of the question. In this sense, control of the population 
transfer branching ratio is incomplete. 

Consider now the null eigenvector |Ai) with Sip / 
and flc 7^ 0. First, if the pump pulse follows the Stokes 
pulse and the control pulse follows the branching pulse, 
then |Ai) will first correlate with the initial state and then 
correlate with the product states as time evolves. Second, 
the two nodes of |Ai) on states |2) and |5) are indepen- 
dent of the phases and amplitudes of the control pulse 
and the pump pulse. Hence, while |Ai) remains the key 
eigenstate for adiabatic population transfer whilst main- 
taining zero population of the intermediate state and the 
branch state, it possesses additional experimentally con- 
trollable parameters that are absent in \0)kr and \0}'kr- 
Third, it is seen that the value of B associated with |Ai) 
can be time-dependent during the population transfer 
and is sensitive to the relative phases of the laser pulses. 
This phase sensitivity is intriguing and makes it clear 
that our five-level four-pulse phase-dependent extended 



The key element in our approach to controlling popu- 
lation transfer is motivated by the observation that the 
laser pulses can be designed such that a third node can 
be introduced to the null eigenvector. From Eq. I|1U|) it 
is seen that if the control fields can ensure that 



O* O* 

1 L pi £^4 



0, 



(14) 



then 



|Ai) = 



B4 



o* 

"_B4 




-n* c 




(15) 



That is, |Ai) will have an additional node on state |3). A 



similar condition i.e., fi^fi? 



O* O* 



will induce a 



node on state |4). With this third node of |Ai) transfer 
of population to one of the degenerate product states can 
be completely shut off. Significantly, since the third node 
on state |3) (if obtained) implies that the time-evolving 
adiabatic state |Ai) will have zero overlap with state |3) 
at all times, it becomes possible that state |3) is never 
populated during the population transfer from the initial 
state to the target state |4). 

Condition l|14l) can be met at all times by controlling 
the phases and the envelope functions of the four laser 
pulses. As the Rabi frequencies are complex in general 
and time-dependent, condition ljl4|l requires matches of 
both the laser phases <px, and their envelope functions 
fx(t). Specifically, condition (fT4*jl holds if 



+ 4>c - 4>p — 4>s 



«45 



Q:24 + «12 — C*15 
s |2~ 



E P E B \^i2\\Hi5 \ exp 



(tp - t B y 



T 2 



E S E C |^24 1 |Mis I exp 



n 



T 2 



T 2 



rp2rp2 
1 P 1 B 



T 2 S T 2 



Tltp 



Tpt B 



Tltc 



T 2 c t s 



n 



n 



(16) 



(17) 



(18) 



(19) 
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Clearly, with the intensity and phases of all other laser 
pulses fixed, Eqs. ((TCI) and ((T7|l can be fulfilled by scan- 
ning the phase and intensity of the control pulse. Equa- 
tions (|18fl and l|19l) are independent of the molecular 
properties and can also be easily satisfied. We now con- 
sider two sample solutions to Eqs. (|18J) and Ijl9|l that will 
be used in the next section. First, suppose that all the 
pulse widths are the same i.e., Tp — T$ = Tb = Tc = T. 
Then Eq. JTBJ is true; and if t s = 0, t P = T, t B = 0, 
then from Eq. 119fl one obtains tc — T. Second, con- 
sider another example where the branching pulse sits in 
the middle of the pump and the Stokes pulses i.e., ts = 0, 
t P = T,t B = 0.5T, and T B = V2T, T P = T s = T. Then 
one obtains T c = V2T from Eq. JXEJ and tc = 2.5T 
from Eq. I|19|) . Note that these sample solutions have 
ensured that the pulse order is counter-intuitive so that 
|Ai) initially has unity overlap with state |1). Detailed 
conditions for completely suppressing population trans- 
fer to state 1 4) can also be found by slightly modifying 
Eqs. OH and l(T7jl, 

There is only one restriction imposed on the transition 
dipole moments in order to successfully apply the five- 
level four-pulse approach outlined above. In obtaining 
Eq. I|10[l we have assumed that £ls4^l*B3 7^ ^S3^*ba- K 
ils4^*B3 = ^S3^*B4 tnen E( l- 09 suggests that |Ai) has 
a node on state In this case |Ai) cannot be used 
for adiabatic population transfer since it cannot corre- 
late with the initial state. Indeed, if we now further im- 
pose condition l|14f> . then one obtains a third relation 
n* p £l B3 — fl S 3ft c . According to Eq. lfTT|) . this acciden- 
tal case then results in two additional null eigenstates 
that will spoil the branching ratio control. To avoid 
^S4^B3 = ^S3^B4> one obtains from Eqs. ©-© the 
restriction 

1/^23 1 ,., si , l/J-35 1 r xi , on x 

f exp[z(a 23 - a.24)\ T 1 r exp[-i(a 35 - a 45 )J- (20) 

|A*24| |M45| 

Note that the above inequality also guarantees that if 
| Ai ) has a node on state |3) (or |4)) it will not have a 
node on state |4) (or |3)). 

It should be stressed that given Eqs. (|16fl - l|20[l none 
of the three states |2), |3), and |5) will, in the adia- 
batic limit, be populated during the population transfer, 
and yet there is 100% population transfer from the ini- 
tial state to state |4). Consider now more realistic cases 
where the lifetimes of the product states are short com- 
pared with the duration of the laser pulses. We assume 
that the control fields are sufficiently strong that the as- 
sociated peak Rabi frequencies are large compared with 
the decay rate constants of the product states. Then 
the lifetime effects can be described by a first-order per- 
turbative treatment [l2l l30| | . In the adiabatic limit this 
gives 

P 3 = 2r 3 ^ t/ |<A 1 |3)| 2 dt, (21) 

where T 3 is the decay rate constant of state |3) and P 3 
is the associated yield. Equation l|21|) indicates that if 



the adiabatic state |Ai) populates state |3) at some in- 
termediate time, then P 3 is necessarily nonzero. That 
is, if state |3) decays, then the mere condition that the 
adiabatic state |Ai) does not overlap with state |3) at 
later times does not suffice to achieve complete suppres- 
sion of P 3 . However, if |Ai) has a node on state |3), 
then (Ai|3) = at all times and Eq. lj?T)) gives P 3 = 0. 
Clearly, then, condition i|14|) for the complete suppres- 
sion of P 3 is also valid when the product states decay, 
and their decay rate constants are relatively small com- 
pared with the various peak Rabi frequencies (for exam- 
ple, if the lifetime of the target states is of the order of 
100 femtoseconds, then the required peak laser intensity 
should be of the order of 10 11 W/cm 2 ). That is, since 
the unwanted state |3) can be never populated, then the 
complete suppression of its yield can also be achieved 
even when the lifetimes of the product states are short 
compared with the duration of the laser pulses. 

The level configuration shown in Fig. ^> resembles 
that in the well-known "two-photon + two-photon" co- 
herent phase control scenario |2l l3S| . However, the asso- 
ciated control mechanisms are drastically different. In 
accord with the weak-field physical picture associated 
with the 'two-photon + two-photon" approach, one intu- 
itively expects that if |Op| • |^5 3 | = \$lc\ ■ |^b 3 | then 
the interference between the two excitation pathways 
|1) -> |2) -> |3) and |1) -► |5) -> |3) is most effective 
and the phase control will be optimal. This intuition 
is irrelevant and useless here, since condition l|14|) is a 
remarkable strong-field result and is extremely counter- 
intuitive. For example, condition l|14f) . which is for the 
complete suppression of the yield of state |3), does not 
specifically contain any information about the transition 
dipole moments that are related to state |3). In this 
sense the control is universal: as long as conditions 114|) 
and H20|) are met and the dynamics is adiabatic, a state 
that is degenerate with, and orthogonal to, state |4) is 
guaranteed to have a zero yield, even if detailed infor- 
mation about the properties of this unwanted state is 
unavailable. Another important difference between the 
"two-photon + two-photon" scenario and our five-level 
four-pulse scheme is that the former can be simplified in 
two special cases. Specifically, if the energy difference 
between states |5) and |1) is identical with that between 
states 1 4) and |2), or if state |5) is degenerate with state 
|2), then only two laser pulses are needed to realize the 
weak-field "two-photon + two-photon" coherent control. 
However, using only two laser pulses in these two special 
cases is not feasible here due to the conditions required 
for adiabatic passage [see Eqs. (|T5|l . (|TT))l - ijlT?)l ]. 



III. CREATION OF ARBITRARY 
SUPERPOSITION OF TWO DEGENERATE 
TARGET STATES 

Of course, states |3) and |4) are just two candidate 
basis states of a two-dimensional degenerate subspace. 
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Consider now two new and arbitrary orthogonal basis 
states of this subspace, 

|3') = sin(0)|3) + exp(i/3)cos(0)|4), (22) 
|4') = cos(0)|3) -exp(i/3)sin(0)|4), (23) 

where 9 and (3 are two independent parameters character- 
izing the superposition states |3') and |4'). The transition 
dipole moments related to states |3') and |4'), denoted by 
M231 M24j A*35j an d can be obtained from linear com- 
binations of ^23 j M24, M35J and M45- Specifically, 



M23 



M24 



a4 



M45 



sin(0)exp(«a23)|/X23| 

+ exp[i(/? + a 24 )]cos(0)| M24 |, (24) 

cos(0)exp(ia 23 )|^ 23 | 

-exp[i(/? + a 24 )]sin(0)|/i 24 |, (25) 

sin(0)exp(ia 35 )|^35| 

+ exp[-i(/3-a45)]cos(0)|^45|) (26) 
cos(0)exp(za 35 )|^35| 
-expH(/3-Q!45)]sm(0)| M 45|. (27) 



Now if we apply the formalism outlined in the previous 
section to states |3') and |4') instead of states |3) and 
|4), we find that the yield of state |3') can be completely 
suppressed, resulting in 100% population transfer to state 
|4'). Not surprisingly, the specific conditions for the com- 
plete suppression of the yield of state |3') are the same as 
Eqs. O-lJTSJ, except that |^ 45 | — > |/4 5 |, IM24I — > IM24L 
Q!24 — > «24i an< i a 45 ~~ * a 45i where 

| M24 | - v/cos^^l^p +si n 2 (0)| M24 |2, (28) 



I /4s I 



cos(0) sin(a 2 3)|^ 23 | 



cos(0) cos(a 23 )|^ 23 | 
cos(0) sin(a 35 )|/i3 5 | 



2 (0)|^ 35 | 2 +sin 2 (0)|M45| 2 , 
tan(a 24 ) = 
- sin(0) sin(/? + a 24 )|^ 24 | 



- sin(0) cos(/3 + a 24 )|/i 24 | ' 
tan(a 45 ) = 
sin(0) sin(^ - a 45 )|^ 45 | 



cos(0) cos(a35)|/i 35 | - sin(0) cos(/3 - a 45 )|^ 45 | 



(29) 



(30) 



(31) 



Thus, an arbitrary superposition of two degenerate target 
states, such as state |4') in Eq. can be created with 

our five-level four-pulse scheme if the lifetimes of states 
|3) and |4) are negligible compared with the duration 
of the laser pulses. Evidently, this also indicates that 
an arbitrary population transfer branching ratio between 
two degenerate product states can be attained. 

Our approach to the creation of an arbitrary superpo- 
sition of two degenerate target states has two important 
advantages. First, the unwanted superposition state |3'), 
as discussed above, can never be populated because the 
control is realized via the time-evolving adiabatic state 
that has a node on state |3'). Second, it is assumed that 
any single laser field cannot individually address the de- 
generate target states. In contrast, previous proposals 
flH fTsL \yA [H| for the creation of superposition states 



of magnetic sublevels have at most only one of the two 
features. 

The ability to create an arbitrary superposition of two 
degenerate molecular states is of considerable theoreti- 
cal and experimental interest. Suppose state |3) is as- 
sociated with one potential energy surface and state |4) 
is associated with another potential energy surface, i.e., 
|3) = |n 3 ) ® |e 3 ) and |4) = |n 4 ) ® |e 4 ), where |n 3 ) and 
|n 4 ) describe the rovibrational motion of states |3) and 
|4), and |ea) and |e 4 ) denote the associated electronic 
states. One then obtains that the superposition state 
|4') is an entangled state: 

|4') = cos(0)|n 3 ) ® |e 3 ) - sin(0) exp(i/3)|n 4 ) ® |e 4 ). (32) 

As such, via the creation of an arbitrary superposition 
of two degenerate target states, our control scheme also 
offers a means for the creation of arbitrary quantum en- 
tanglement between the rovibrational and electronic de- 
grees of freedom. The arbitrary quantum entanglement 
in molecular systems thus created may be of importance 
for quantum information science 39]. In addition, if the 
created superposition state |4') can be used as an initial 
state in other processes, such as for molecular scatter- 
ing, then the two participating degenerate states |3) and 
|4) will be able to interfere with each other. The asso- 
ciated quantum interference effects can be manipulated 
by changing the parameters 9 and /?, giving rise to an 
opportunit y fo r quantum interference control of collision 
events [ioll^. 



IV. NUMERICAL EXAMPLES 

In this section we provide some numerical examples of 
the control of population transfer in our five-level four- 
pulse scheme. Consider first the results shown in Figs. 
2 and 3 where the control goal is to completely suppress 
the yield of state |3). To obtain the results shown in Fig. 
2, we assumed that the product states do not decay and 
that 

n P T = 40/p(t), (33) 
n Si T = {20 + 20i)f s (t), (34) 
n Bi T = (30 + 200/sW, (35) 

n c T = (10 + 50i)/c(*), (36) 



with 



fp{t) = fc(t)=exp[-(t-T) 2 /T% (37) 
fs(t) = f B (t)=eM-t 2 /T 2 ). (38) 

With these choices it is seen that condition l|14|) is met, 
so the yield of state |4) is expected to be almost 100% 
and the yield of state |3) is expected to be almost zero, 
irrespective of the properties of state |3). Two cases 
with drastically different transition dipole moments that 
are related to state |3) are considered. In case (a) 
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FIG. 2: Numerical examples of the dynamics of adiabatic 
population transfer in the five-level four-pulse phase-sensitive 
extended STIRAP scheme. Shown here is the time depen- 
dence of the population on state \k) (denoted P&). QpT = 
40/p(t), n S 4T = (20 + 20i)f s (t), VL Bi T = (30 + 20i)f B (t), 
QcT = (10 + 50t)/ o (t), with f P (t) = f c {t) = exp[-(t - 
T) 2 /T 2 ] and f s (t) = f B (t) = exp(-i 2 /T 2 ). In case (a) 
QssT = 3Q/ 5 (*) and tt m T = 20/s(t); in case (b) fi S3 T = 
10/s(t) and Q B3 T = (80i)f B (t). Note that during the popu- 
lation transfer P%, P3, and P5 are negligible at all times. 



transfer will be further decreased by more than one order 
of magnitude. In addition, for given field intensity and 
pulse duration one may optimize the time delay between 
the laser pulses to obtain the smallest nonadiabaticity. 
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FIG. 3: Same as in Fig. 2 except that flcT and the 
pulse envelope functions are different. Here £lcT — (10 + 
50i)/c(t)exp(2.0), f P (t) = exp[-(i - T) 2 /T 2 }, fs(t) = 
exp(-t 2 /r 2 ), f B (t) = exp[-0.5(t - 0.5) 2 /r 2 ], and f c (t) = 
exp[-0.5(t - 2.5T) 2 /T 2 ]. As in Fig. 2, cases (a) and (b) 
assume different properties of state |3). 



il S3 T = 30f s (t) and n B3 T = 20/ B (t), and in case (b) 
il S3 T = 10f s (t) and fl B3 T = (8(H)/ B (i). As seen from 
the detailed time dependence of the population (denoted 
Pk) on state \k) (k = 2 — 5) shown in Fig. 2, the results 
in both cases (a) and (b) confirm our previous analy- 
sis. In particular, the population transfer to state |4) is 
complete, the maximal P 3 during the population trans- 
fer is less than 0.05% in case (a) and less than 0.07% 
in case (b), and in both cases Pi and P5 are less than 
0.4% at all times. Of course, due to some weak nonadia- 
batic effects associated with the dynamics of population 
transfer, P2, P 3 , and P5 are not perfectly zero. As in the 
three-level version of STIRAP, using stronger laser fields 
or larger pulse durations will greatly enhance adiabatic- 
ity and therefore further suppress P2, P3, and P5. For 
example, we have checked that if all the Rabi frequen- 
cies or the pulse widths are increased by a factor of ten, 
then the maximal P2, P3, and P5 during the population 



Figure 3 displays the analogous results for a different 
choice of the pulse envelope functions. Here the branch- 
ing pulse is in the middle of the pump and the Stokes 
pulses, and the control pulse is the most delayed pulse. 
As discussed above, condition (fHIJl can also be easily met. 
Except for Vt c T = (10 + 50i) f c (t) exp(2.0), other Rabi 
frequencies here are chosen to have the same peak values 
as those used in Fig. 2, with their envelope functions 
given by 

f P (t) = cxp[-(i- Tf/T 2 ], (39) 
f s (t) = ex P H 2 /T 2 ), (40) 
f B (t) = cxp[-0.5(t-0.5T)7T 2 ], (41) 
fc(t) = exp[-0.5(t-2.5T)7T 2 ]. (42) 

As seen in Fig. 3, the control obtained is almost perfect, 
for both cases (a) and (b) that assume different properties 
of state |3). 
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FIG. 4: Same as in Fig. 2 except that the degenerate product 
states 1 3} and |4) have a common finite lifetime that is ten 
times smaller than the pulse duration. 
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FIG. 5: Dynamics of adiabatic population tranfer from the 
initial state |1) to a superposition of two degenerate target 
states given by |4'} = (|3) + i|4))/V2. P41 and Py represent 
the population on states |4'} and |3'} = (|3) - i\4))/y/2, P2 
and P5 represent the population on states |2) and |5). The 
pulse envelope functions are given by Eqs. tSTl and <I351. 
flc T = (80 + A0i)fc(t), and all other Rabi frequencies are 
chosen to be the same as those used in case (a) in Fig. 2. 



Given these conditions, we find that in order to induce 
a node on state |3') the control pulse should be designed 
to give fl c T = (80 + A0i)f c (t). As shown in Fig. 5, this 
indeed gives almost 100% yield of the superposition state 
1 4'), with the population on states |2), |5), and |3') being 
negligible at all times. Similar results are obtained for 
other superposition states as target states. 



As mentioned above, the almost complete suppression 
of P3 at all times suggests that the control of population 
transfer can hold when the product states rapidly decay. 
This is also confirmed numerically. For example, in Fig. 
4 we show the population dynamics when both the prod- 
uct states have a lifetime that is one order of magnitude 
smaller than the pulse duration, with all other parame- 
ters same as in Fig. 2. It is seen that although the final 
value of P3 is larger than in the non-decaying case in Fig. 
2, it is still less than 0.3%. It is also seen that there is an 
increase in P2 and P5 at intermediate times (but they are 
still negligible). This result is understandable because in 
the presence of decaying product states the dynamics is 
not fully coherent. The situation can always be further 
improved by use of stronger laser fields. 

We have also numerically examined our approach to 
the creation of an arbitrary superposition of states |3) 
and 1 4) whose lifetimes are much longer than the pulse 
duration. One computational example is shown in Fig. 
5, where the target state is |4') = (|3) + i\4))/s/2. Then 
the unwanted superposition state is |3') = (|3)— 1|4))/\/2. 
The pulse envelope functions are given by Eqs. (|37(l and 
iJSSJ, and except for fi^, all other Rabi frequencies are 
chosen to be the same as those used in case (a) in Fig. 2. 



V. CONCLUDING REMARKS 

Using a five-level model, we have shown that complete 
control of the population transfer branching ratio be- 
tween two degenerate product states is achievable using 
four laser pulses with carefully designed relative phases 
and envelope functions. The complete control of pop- 
ulation transfer stems from a peculiar light-matter in- 
teraction in the presence of strong laser fields: the con- 
trol mechanism lies in a three-node null eigenstate that 
can correlate with an arbitrary superposition of two de- 
generate target states. Potential applications of our re- 
sults include absolute selectivity of product formation in 
a photochemical reaction and the creation of an arbitrary 
superposition of two degenerate target states. 

The most fascinating and somewhat surprising feature 
of our five-level four-pulse scheme is that it allows for 
the prediction of complete suppression of the yield of one 
unwanted state of two degenerate product states without 
studying its properties beforehand. This feature of the 
scheme likely has some deep implications. On the other 
hand, it becomes very interesting to ask if there exists a 
control scheme that can guarantee 100% yield of one of 
two degenerate product states with its main properties 
(such as its transition dipole moments) unknown to us. 
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We are not yet ready to give a conjecture here, but believe 
that if the answer is no then the reason could be related 
to the quantum no-cloning theorem 

This study also suggests that STIRAP-like dynamics 
can shed considerable light on, and may provide a sys- 
tematic solution to, the formal question of controllability 
in a multi-level quantum system without or with con- 
straints. In particular, the availability of a complete pop- 
ulation transfer pathway that never populates a certain 
number of states can be examined by studying the ex- 
istence of a single multi-node dressed eigenstate of the 
system in the presence of control fields. If, by designing 
the control fields, such a multi-node dressed eigenstate 
can be made to exist and can correctly correlate with the 
initial and target states, then the controlled dynamics in 



the adiabatic limit gives a specific physical solution to 
the generation of complete population transfer. 

Previous quantum control approaches based on 
STIRAP-like dynamics are independent of laser phases. 
As shown in this work and in Ref. 01, this is not al- 
ways the case. Rather, one may take advantage of the 
laser-phase-dependent quantum interference between dif- 
ferent STIRAP-like sub-processes to achieve more dra- 
matic goals in the laser control of atomic and molec- 
ular processes. It is our hope that in the near future 
extensions of our five-level four-pulse extended STIRAP 
scheme can provide us other amazing opportunities for 
quantum control. 

This work was supported by the National Science 
Foundation. 
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